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Abstract. The 4.62^tm (2164.5 cm -1 ) 'XCN' band has been detected in the M-band spectra of 34 deeply embedded 
young stellar objects (YSO's), observed with high signal-to-noise and high spectral resolution with the VLT-ISAAC 
spectrometer, providing the first opportunity to study the solid OCN - abundance toward a large number of low- 
mass YSO's. It is shown unequivocally that at least two components, centred at 2165.7 cm -1 (FWHM = 26 cm -1 ) 
and 2175.4cm" 1 (FWHM = 15cm" 1 ), underlie the XCN band. Only the 2165.7-component can be ascribed to 
OCN - , embedded in a strongly hydrogen-bonding, and possibly thermally annealed, ice environment based on 
laboratory OCN - spectra. In order to correct for the contribution of the 2175.4-component to the XCN band, a 
phenomenological decomposition into the 2165.7- and the 2175.4-components is used to fit the full band profile 
and derive the OCN - abundance for each line-of-sight. The same analysis is performed for 5 high-mass YSO's 
taken from the ISO-SWS data archive. Inferred OCN - abundances are < 0.85 % toward low-mass YSO's and 
< 1 % toward high-mass YSO's, except for W33 A. Abundances are found to vary by at least a factor of 10-20 and 
large source-to-source abundance variations are observed within the same star-forming cloud complex on scales 
down to 400 AU, indicating that the OCN - formation mechanism is sensitive to local conditions. The inferred 
abundances allow quantitatively for photochemical formation of OCN - , but the large abundance variations are 
not easily explained in this scenario unless local radiation sources or special geometries are invoked. Surface 
chemistry should therefore be considered as an alternative formation mechanism. 

Key words. Astrochemistry - line: identification - line: profile - molecular data- methods: data analysis - ISM: abun- 
dances -ISM: lines and bands - infrared:ISM 



1. Introduction 

The 4.62 /mi (2165 cm" 1 ) feature, commonly referred to 
as the XCN band, was fi rst detected towar d the mas- 
sive protostar W33A by ISoifer et all (|l979j) . Its pres- 
ence was the first observational indication that complex 
chemistry could be occurring in interstellar ice mantles. 
Subsequently, a similar feature was observed around a 
number of other sources, mostly high-mass young stel- 



lar objects (YSO's') llTeeler et al.ll993T.ll995t iDemvk et al 



1998:IPendleton et al.ll999UOibb et alT2000 



2001 



Whittet et all 1200 ljh one field star llTegler et al 



Keane et al 



1995) and several gala ctic centre sources IjChiar et al 
20021 ISpoon et alJl200: " 
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with the participation of ISAS and NASA. 



IPontoppidan et alJ (|2003^ (henceforth referred to as 
Paper I) have recently presented M-band spectra of 44 
YSO's, 31 of which are low-mass sources (<2M©, < 
50L Q ). Most of these sources are deeply embedded in dark 
clouds, with few external sources of ultraviolet radiation, 
offering a unique opportunity to study the chemical char- 
acteristics of ice mantles in low-mass star-forming envi- 
ronments. In addition to the strong CO-ice band, observed 
on most lines of sight, Paper I reported the detection of a 
weaker band in 34 of the spectra (27 of which were low- 
mass objects). This band was labelled the 2175 cm -1 band 
because its peak-centre position ranged from 2162-2194 
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cm -1 , with an average value of 2175 cm . These were 
the first reported detections toward a large number of low- 
mass YSO's of a feature similar to the XCN band. 

The XCN band has been studied extensively in 
the laboratory where it is easily reproduced by pro- 
ton irradiatio n llMoore et al.1 Il983h. vacuum ultravio- 
let photolysis llLacv et al.lll984UGrim fc Greenberdll987l 
[ Demyk et al .1 1 19 9 8ft . or therm a l anne aling l|Raunieret*al] 
120031 Ivan Broekhuizen et alJ l2004h of interstellar ice 
analogues. Experimental studies using isotopic sub- 
stitution proved unequivocally that the carrier of 
the l a boratory feature is OCN~ llSchutte k. Greenberd 
1997t iBernstein et all l200(j iNovozamskv et alJ l200lT 
Palumbo et al.ll2000|) . Despite this assigment, alternative 
identifications of the in terstellar feature are still debated 
l|Pendleton et al.lll999j) . 

Here we present a detailed analysis of the 2175 cm -1 
feature (henceforth described in this paper as the XCN 
band) for 34 of the YSO's from Paper I, plus 5 addi- 
tional high-mass YSO's taken from the ISO-SWS data 
archive, i.e. AFGL2136, NGC 7538 IRS1, RAFGL 7009S, 
AFGL989 and W33A. Our analysis suggests that two 
components are required to reproduce the observational 
line profiles of the XCN band, only one of which can be 
associated with OCN~. The carrier of the second com- 
ponent is not known but may be d ue to CO gas-grai n 
interactions, as described in detail bv lFraser et alJ l|2005h . 

In Sect. 2, the observational details are summarised, 
followed by an overview, in Sect. 3, of laboratory spec- 
troscopy of OCN~ in interstellar ice analogues and studies 
of OCN - formation mechanisms under interstellar con- 
ditions, focusing on UV- and thermally induced OCN - 
formation. The fitting methods used for the XCN band 
analysis are described in Sect. 4, and the results presented 
in Sect. 5. In Sect. 6 OCN - abundances are determined 
and the implications for the OCN - formation mechanisms 
toward the low-mass YSO's studied here are discussed. 



2. Observations 

2.1. Observational details 

This work uses the M-band spectra of a large sam- 
ple of deeply embedded young stars, which were first 
presented in Paper I. The data are part of a 3-5 /im 
spectroscopic survey of low-mass embedded objects in 
the nearest star-forming clouds (p Oph, Serpens, Orion, 
Corona Australis, Chamaeleon, Vela, and Taurus) using 
the Infrared Spectrometer And Array Camera (ISAAC) 
mounted on UTl-Antu of t he Very Large Telescope (VLT) 
( van D ishoeck et al . 2003). All sources have spectral en- 
ergy distributions representative of Class I objects, with 
typical lifetimes of ~ 10 5 yr since cloud collapse. The M- 
band spectra were obtained in the 4.53-4.90 /im (2208- 
2040 cm -1 ) range using the medium resolution mode, re- 
sulting in a resolving power of A/AA = 5000-10000. The 
reduction of the spectra is described in detail in Paper 
I. In short, the spectra were corrected for telluric ab- 
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Fig. 1. (a) Observed peak-centre positions of the XCN 
band toward all 39 low- and high-mass YSO's studied 
here. Bands originating in high-mass YSO spectra are rep- 
resented in black, all other YSO's (dominated by low-mass 
objects) are hatched. Peak centres are binned at 1.6 cm -1 . 
(b) Peak-centre positions of laboratory spectra of the 1/3 
vibrational band of solid OCN - in a variety of interstel- 
lar ice analogues, also summarised in Tables ^ and [21 
H 2 0- and CH 3 OH-rich ices are shown in black, H 2 0- and 
CHaOH-free ices are hatched, and those thermally an- 
nealed to > 110 K are white. These centre positions are 
rounded to the nearest integer wavenumber (cm -1 ). 



sorption, flux-calibrated relative to bright early-type stan- 
dard stars, and wavelength-calibrated relative to telluric 
absorption lines in the standard star spectrum with an 
accuracy of ~ 5 km s . The final optical depth spec- 
tra were derived by fitting a blackbody continuum to the 
4.52-4.55/nn (2212-2198cm -1 ) and 4.76-4.80/im (2101- 
2083 cm -1 ) regions, where no features are expected, tak- 
ing care to exclude known gas-phase lines from the fit. The 
short wavelength end of the spectra is more noisy due to 
the onset of the strong telluric CO2 features. 
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Fig. 2. The distribution of centre positions of the XCN 
band toward the various lines-of-sight in bins of 1.6 cm -1 
including the uncertainties derived from the Monte-Carlo 
simulations. The Gaussian at 2168.7cm -1 is obtained 
from a single Gaussian fit to the full data set, that at 
2166.9 cm -1 from a similar fit excluding all XCN band 
centres > 2172 cm -1 . This analysis strengthens the conclu- 
sion that two components are present in the observational 
data. The dashed line indicates the maximum wavenum- 
ber beyond which no laboratory-based f3(OCN") feature 
has been found in interstellar ice analogues. 



Fig. 3. Example of a laboratory spectrum of v% (OCN ), 
obtained after thermal annealing to 120 K of an ice com- 
posed of H 2 0/HNCO/NH 3 = 120/1/10 (grey curve). 
After accounting for the baseline (dotted line), an excel- 
lent fit to the laboratory data is obtained with a single 
Gaussian fit, centred at 2165.9cm- 1 , FWHM = 23cm- 1 
(black curve). 



In addition, spectra of 5 high-mass sources taken from 
the ISO-SWS data archive, i.e. AFGL2136, NGC7538 
IRS1, RAFGL 7009S, AFGL989 and W33A, were added 
to the sample. 

2.2. The XCN band 

In addition to the CO-ice band, a weak secondary feature 
was detected in a significant subset of the lines of sight. In 
Paper I this feature was reasonably well reproduced using 
a single Gaussian whose peak-centre position varied sig- 
nificantly (from 2162 to 2194cm- 1 ) and Full Width Half 
Maxima (FWHM) ranged from 9 to 36 cm- 1 . This distri- 
bution is summarised in Fig. where the peak-centres 
of all XCN bands in the observed sample are plotted, to- 
gether with the data from the ISO sources. 

Careful inspection of the XCN band shows that in 
some sources like W33 A the feature peaks at 2165.7cm- 1 
whereas in other sources like Elias 32 and IRS 63 the fea- 
ture is shifted much further to the blue, around 2175 cm- 1 . 
This conclusion that two components underlie the XCN 
band is supported by a statistical analysis of the observed 
centre positions shown in Fig. [21 The uncertainties in the 
observed distribution were quantified using a Monte Carlo 
approach. Each Monte Carlo run varies the XCN band 
centre position of each source within a normal distribu- 
tion with a standard deviation from peak-centre position 
as found in Paper I. This produces a perturbed distribu- 
tion allowed within the uncertainties on the band posi- 
tions. By calculating a large number (10 5 ) of perturbed 
distributions, an uncertainty on the number of sources 
within each wavenumber bin can be estimated by sim- 



ply deriving the standard deviation for the set of per- 
turbed distributions. A single Gaussian distribution cen- 
tred at 2168. 7± 1.9 cm- 1 (x 2 = 5.58) is seen to give a 
poor fit. When all bands centred at > 2172 cm- 1 are ex- 
cluded from the fit, the centre of the Gaussian distribution 
shifts to 2166.9 ± 1.5cm- 1 (x 2 — 1.57). The converging x 2 
confirms that two components, i.e. one centred at around 
2166cm- 1 and one at around 2175cm- 1 , provide a much 
better fit. 



2.3. fyO-ice column densities 

In order to derive the OCN- abundances, water ice col- 
umn densities were estimated using L-band spectra (2.85- 
4.2 /im) from the same survey (Dartois et al., in prep.). 
Optical depth spectra of the H^O-ice feature were ob- 
tained by fitting a blackbody spectrum to a -fT-band pho- 
tometric point from the 2MASS point source catalogue 
and the available spectroscopic points in the L- and M- 
band spectra at wavelengths longer than 3.8 /im. Care was 
taken to exclude any hydrogen recombination lines from 
the fit. Estimated uncertainties in the water ice optical 
depths are ~ 20%, mainly due to the vari able near-in frared 
fluxes of embedded young stars (see e.g. iKaasIl 19991) . The 
optical depths were converted to column densities using a 
factor of 1.56 x 10 18 cm -2 . This factor was derived by in- 
tegrating over a laboratory-based H2O spectrum fitted to 
the highest quality water bands a nd scaling with a ban d 
strength of 2 x lO -16 cmmolec. -1 l|Gerakines et, aljfl995h . 
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Table 1. Laboratory experiments on the spectroscopy of ^(OCN - ) in interstellar ice analogues (T = 10-100 K) a 



T 

(K) 



10 

10 

10 

10 

10 

10 

12 

IS 

IS 

IS 

IS 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

SO 

30 

30 

35 

50 

50 

55 

60 

80 

80 

80 

80 

100 

100 

100 

100 

100 



Energetic 
processing 



Ice composition 
(spectral resolution where known /cm" 1 ) 



^ 3 (OCN-) 



ion (30kcV) 
thermal 
UV (at 10K) 
UV (at 10K) 
UV (at 10K) 
UV (at 10K) 
UV (at 12K) 
UV (at 12K) 
UV (at 10K) 
UV (at 3 OK) 
ion (60keV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
proton (0.8MeV) 
UV (at 15K) 
UV (at 15K) 
UV (at 1SK) 
UV (at 15K) 
proton (lMcV) 
proton (lMcV) 
proton (lMcV) 
proton (lMeV) 
proton (lMcV) 
proton (lMcV) 
proton (lMcV) 
proton (lMeV) 
proton (lMcV) 
proton (lMeV) 
proton (lMcV) 
proton (lMcV) 
proton (lMcV) 
UV (at 12K) 
thermal 
UV (at 10K) 
proton (0.8 MeV) 
UV (at 12K) 
UV (at 37K) 
UV (at 10K) 
UV (at 12K) 
UV (at 12K) 
UV (at 12K) 
UV (at 12K) 
UV (at 12K) 
UV (at S7K) 
UV (at 37K) 
ion ( 60ke V) 
UV (at 10K) 
UV (at 10K) 



H 2 0/CH 4 /NH 3 = 5:4:2 (2) 
HNCO/NH 3 = 1:10 (0.5) 
H 2 0/CO/CH 4 /NH 3 = 6:2:1:1 (2) 
CO/NH 3 = 3:1 

H 2 0/CO/CH 3 OH/NH 3 = 100:10:50:10 

H 2 0/CO/CH 3 OH/C„H 2 „ +2 /NH 3 = 100:10:50:10:10 (0.9) 
HNCO/NH 3 = 1:1 (2) 
HNCO/NH 3 = 1:100 (2) 
CO/NH 3 = 1:40 (2) 
CO/NH 3 = 1:40 (2) 
C0 2 /N 2 = 1:1 (1) 
H 2 0/CO/NH 3 = 5:1:1 (4) 
H 2 0/C0 2 /NH 3 = 1:1:2 (4) 
H 2 0/CH 4 /NH 3 = 1:1:1 (4) 
H 2 0/C0 2 /N 2 = 5:1:1 (4) 
H 2 0/CH 4 /N 2 = 1:1:1 (4) 
CO/NH 3 = 1:1 (4) 
C0 2 /NH 3 = 1:2 (4) 
CO/N 2 = 1:1 (4) 
CO/N 2 = 1:2 (4) 
H 2 0/HNCO/NH 3 = 120:1:10 (2) 
H 2 0/HNCO/NH 3 = 120:1:10 (2) 
HNCO/NH 3 = 1:10 (2) 
H 2 0/HNCO/NH 3 = 140:8:1 (2) 
H 2 0/CH 4 /NH 3 = 1:2:3 
H 2 0/CH 4 /NH 3 = 1:2:4 
H 2 0/CH 4 /NH 3 = 2:4:7 
H 2 0/CH 4 /NH 3 = 1:7:10 
H 2 0/CH 4 /NH 3 = 1:3:6 
H 2 0/CH 4 /NH 3 = 1:3:3 
H 2 0/CH 4 /NH 3 = 1:4:3 
H 2 0/CH 4 /NH 3 = 15:7:9 
H 2 0/CH 4 /N 2 = 1:1:1 
H 2 0/CO/N 2 = 5:1:1 
H 2 0/CO/NH 3 = 2:2:1 
H 2 0/CO/NH 3 = 2:1:3 
H 2 0/CO/NH 3 = 5:1:10 
HNCO/NH 3 = 1:100 
HNCO/NH 3 /Ar = 2:2:1000 (0.125) 



H 2 0/CO/CH 3 OH/C„H 



/NH 3 = 100:10:50:10:10 (0.9) 



CO/CH 4 /N 2 = 1:1:100 (1) 
HNCO/NH 3 = 1:100 (2) 
CO/NH 3 = 1:1 (2) 

H 2 0/CO/CH 3 OH/C„H 2 „ +2 /NH 3 = 100:10:50:10:10 (0.9) 

HNCO/NH 3 = 1:100 (2) 

CO/NH 3 = 1:1 (4) 

H 2 0/CO/NH 3 = 1:1:1 (4) 

H 2 0/CO/NH 3 = 2:1:1 (4) 

H 2 0/CO/NH 3 = 3:1:1 (4) 

CO/NH 3 =1:1 

CO/NH 3 = 1:1 (2) 

C0 2 /N 2 = 1:1 (1) 

H 2 0/CO/CH 3 OH/NH 3 = 100:10:50:10 

H 2 0/CO/CH 3 OH/C, 1 H 2 „ +2 /NH 3 = 100:10:50:10:10 (0.9) 



cm -1 ) 


(cm" 1 ) 




2170 




1 


2151 




2 


2167 


25 


3 


2155 




4 


2168 




5 


2160 




6 


2155 




7 


2155 




7 


2152 




7 


2152 




7 


2168 




8 


2167 


26 


9 


2164 


25 


9 


2159 


23 


9 


2170 


25 


9 


2159 


27 


9 


2164 


27 


9 


2158 


26 


9 


2164 


27 


10 


2158 


26 


10 


2169 


25 


11 


2171 


25 


11 


2160 


20 


11 


2165 


25 


11 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2170 




12 


2155 




7 


2157 




2 


2163 




6 


2166 


25 


13 


2154 




7 


2158 




7 


2163 




6 


2154 




7 


2157 




14 


2160 




14 


2163 




14 


2166 




14 


2163 




15 


2157 




7 


2165 




8 


2166 




5 


2163 




6 



"Experiments i n italics show ices that lack H2 O. All 1/3 peak-posit ion s are rounded up o r down t o the nearest inte ger 
wavenumber . 1 IStrazzulla fc Palumbol Jl998l). 2 lRaunier et all J2003h 3 Id'Hendecourt et al] lll986ft. 4 lLacv et alJ Jl984h. 5 
Teeler et alJ il993h. 6 lAllamandola et al.l dl988l). 7 lNovozamskv et al.l feoOlft. 8 |Palumbo et al.l feOOCh. 9 iHudson fc Moord 



(2000) 



Hudso n et al] (1200 ll). 11 Iv an Broekhuiz en et al] yOOJ), 12 iMoore et al] J1983T) . 13 iMoore fc Hudson! (2003), 



iGrim fc Greenberd (Il987l) . 15 iBernstein et al] fcOO' 



3. Laboratory experiments on OCN pertinent to 
the XCN band analysis 

3.1. Spectroscopy 

The formation and infrared spectroscopy of solid OCN~ 
has been studied extensively in numerous laboratories. Its 
spectrum is dominated by the CN stretching-vibration 
(the vn, band), which can be well fitted by a single 
Gaussian profile (see for example Fig. [3Jl. In Tables Q 
andean overview is given of all (to the best of our knowl- 
edge) published peak centres and FWHM of the ^3 band 
of OCN - , as measured in the laboratory under condi- 
tions relevant to the interstellar case, e.g. matrix-isolation 
studies have been excluded. In addition, information is 



tabulated on the initial ice composition, ice temperature 
and processing mechanism applied to the ice. With the 
exception of H2O, CH3OH and NH3, most of the other 
ice components have fully desorbed by 120 K under high- 
vacuum laboratory conditions but OCN - remains in the 
solid state, stabilised by a counter io n (possibly NH4 , 
iDemvk et all Il998l FKban et all l2003|) . Hence, Table [l] 
contains those experiments in which the ice matrix is pre- 
dominantly intact, whereas Table El summarises those in 
which most ice constituents have desorbed or are desorb- 
ing. 

Fig. shows the distribution of peak centres of the 
laboratory 1/3 of the OCN - feature highlighting the initial 
ice environment, i.e. H2O- or CHsOH-rich (black), H2O- 
and CHsOH-lacking (hatched) and thermally annealed, 
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Table 2. Laboratory experiments on the spectroscopy of OCN - in interstellar ice analogues (T > 100K) a 



T 

(K) 


Energetic 
processing 


Ice composition 
(spectral resolution where known /cm -1 ) 


1/3 (OCN ) 
(cm" 1 ) 


FWHM 

(cm- 1 ) 


rcf . 


110 


thermal 


H 2 0/HNCO/Ar = 12:1:1000 (0.12) 


2170 




1 


120 


UV (at 10K) 


CO/NH 3 =1:1 


2155 




2 


120 


UV (at 15K) 


H 2 0/HNCO/NH 3 = 120:1:10 (2) 


2166 


22 


3 


130 


thermal 


H 2 0/HNCO = 10:1 (0.5) 


2172 




1 


140 


UV (at 37K) 


CO/NH 3 = 1:1 (2) 


2159 




4 


ISO 


ion (60keV) 


C0 2 /N 2 = 1:1 (1) 


2163 




5 


150 


UV (at 10K) 


CO/NH 3 = 3:1 


21 65 




6 


150 


UV (at 10K) 


H2O/CO/CH3OH/NH3 = 100:10:50:10 


2162 




7 


ISO 


thermal 


HNCO/NH 3 = 1:10 (0.5) 


2165 




1 


180 


UV (at 12K) 


CO/NH 3 = 1:1 (4) 


2157 




8 


ISO 


UV (at 12K) 


CO/NH 3 = 1:10 (4) 


2157 




8 


200 


ion (oOkeV) 


CO2/N2 = 1:1 (1) 


2162 




5 


200 


UV (at 10K) 


H 2 0/CO/CH 4 /NH 3 = 6:2:1:1 (2) 


2167 


25 


9 


200 


UV (at 10K) 


H2O/CO/CH3OH/NH3 = 100:10:50:10 


2164 




7 


200 


UV (at 10K) 


H2O/CO/CH3OH/NH3 = 100:10:50:10 (0.9) 


2165 




10 


200 


UV (at 10K) 


H20/CO/CH30H/C„H 2 „ + 2/NH3 = 100:10:50:10:10 (0.9) 


2164 




10 


200 


UV (at 10K) 


H 2 0/CO/CH 3 OH/C„H2n+2/ NH 3 = 100:10:50:10:10 (0.9) 


2163 




10 



a Experiments in italics show ices that lack H2O. All 1/3 pe a k-pos itio ns are rounded up or down to the nearest integer 
waven u mber. 1 [k aunier et al] i2003t). 2 IS chutte fc Greenberd lll 997f). 3 Ivan Broekhuizen et alJ l|2004lh *jNovozainsk^^^l 
|200lh . 5 IPalumbo et allfeOOfllh 6 | Lacv et, all l|l984l h 7 lTegler et all (jl99.TT l. 8 lGrim fc Greenberd fi987f l 9 ld'Hendecourt et al. 
Jl986li . 10 lAllamandola et alJ ll98Sft 



i.e. those listed in Tabled (T > 110 K, white). It is impor- 
tant to note that this histogram is significantly weighted 
by the parameter space investigated in the laboratory. 
Peak centres vary strongly with ice composition, ranging 
from 2151 to 2172cm" 1 . Within this range, i/ 3 (OCN - ) 
peaks at betwee n 2151-2160 cm -1 in the absence of 
H 2 Q and CHcOH llNovozamskv et all200lllRaunier et all 
l200.il iGrim fc Greenberd Il987t iHudson fc Moorel l2000l) . 
or at between 2160-2172 cm -1 when associated with 
a strongly hydrog en-bonding environment, i.e. H2O- or 
CH 3 OH-rich ices JGrim fc Greenberdll987|) . No trend is 
apparent in the data relating the actual formation mech- 
anism of OCN - to its band position. 

From Fig. is also clear that in warm ices 
(T > 110 K), zy 3 (OCN - ) peaks at similar positions to the 
colder ices. However within a single ice matrix, thermal 
annealing beyond 100 K tends to shift the ^ 3 (OCN - ) peak 
to betwee n 2155 and 2167 cm" 1 when OCN~ is stabilise d 
by NH4 l|Novozamskv et al.ll200lt IPalumbo et al.ll2000|h 
In the exclusive presence of H2O (i.e. for H2O/HNCO and 
H 2 0/HNCO/Ar), however, j/ 3 (OCN - ) must be stabilised 
by water solvation (most plausibly H30 + ) and seems ther- 
mally unaffected, peaking at 2170 to 2172 cm" 1 . 

The spectral range of peak positions observed for labo- 
ratory vz (OCN - ) implies that only one component of the 
observed XCN-bands can be associated with OCN - . No 
laboratory i/ 3 (OCN - ) spectrum in Fig. ^ peaks beyond 
2172cm -1 , whereas a significant fraction of the astronom- 
ical spectra in Fig. QJi peaks around 2175cm -1 . This is 
emphasised by the dashed line in Fig. [21 indicating that 
only those XCN-bands represented by the Gaussian fit to 
the left of the dashed line can be matched to a laboratory 
spectrum of OCN - . The peak centres of these XCN-bands 
can be explained by OCN - residing in a strong hydrogen- 
bonding ice environment that may have been thermally 
annealed up to 200 K. OCN - embedded in an ice that 
lacks a strong hydrogen-bonding network generally peaks 



red-wards of the XCN-bands observed, even after thermal 
annealing, causing this kind of ice to be an unlikely inter- 
stellar OCN - environment. 

3.2. Band strength 

The reported value of the ^ 3 (OCN - ) ban d strength 
AnriM - , varies considerably in the literature. lLacv et all 
(1983) assumed A OCN - <1.0xl0 -17 cmmolec. -1 , based 
on the assumption that the cross section of a typical CN- 
stretching vibration i s approximately 5 t i mes l ess than 
that of a CO-stretch. Id'Hendecourt et al] l)l986l) derived 
A OCN - from photolysis experiments to be 1.8-3.6x 10 -17 
cmmolec. -1 , assumin g a constant carbon budget per ex- 
periment, followed bv lDemvk et al.l l)l998h . who assumed 
a constant oxygen budget to obtain Aqcn- >4.3x10 -17 
cmmolec. -1 . Both studies probably underestimated the 
value of Aqcn- due to the formation of infrared-weak or 
inactive photoproducts, and the neglect of photodesorp- 
tion of volatile species like CO and CO2 (van Broekhuizen 
et al. in prep). 

More recently, Ivan Broekhuizen et al.l l|2004F) deter- 
mined Aqcn- horn the thermal deprotonation of HNCO 
by NH 3 , using the NH 3 to NH4 conversion as reference, 
giving ^4ocn- = 1.3xl0 -16 cmmolec. -1 . In the present 
analysis of the XCN band, this value of Aqcn- nas been 
adopted because its determination is not influenced by any 
photodesorption or the formation of photo-products (both 
known and unknown), which may have affected the pre- 
viously derived values. The adopted band strength does 
not affect any subsequent conclusions in this paper on 
the photo-production of OCN - , since these reaction ef- 
ficiencies were derived using the same value of the band 
strength. However, the adopted Aqcn- value does impact 
on the conclusions drawn if OCN - is thermally formed, 
because significantly more HNCO will be required in the 
solid state if a smaller value of A OCN - were used. 
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Table 3. Components adopted in the fit of the observed 
XCN bands 



Wavelength [^m] 



Fig. 4. The XCN band as observed in the line-of-sight 
toward five YSO's, ordered from top to bottom from 
the most red-centred to the most blue-centred band, i.e. 
(a) W33A, (b) HH46, (c) Reipurth50, (d) Elias32 
and (e) IRS 63. All spectra (black) are continuum sub- 
tracted (small-dashed line) and the total fit is shown 
(grey curve). The two components underlying the XCN 
band, i.e. the 2165.7-component (dashed curve) and 
the 2175.4-component (dotted curve) are shown. The 
three-component fit to the solid CO band at 4.67/im 
(2139 cm" 1 ) is shown only in (b) by the three grey dash- 
dotted curves. The vertical dashed line at 4.62 fim guides 
the eye. 





2165.7-component 


2175.4-component 


Line profile 


Gaussian 


Gaussian 


Centre (cm -1 ) 


2165.7 


2175.4 


FWHM (cm" 1 ) 


26 


15 



3.3. Formation mechanisms 

The presence of OCN - in interstellar ices has become 
widely cited as a probe of energetic processing of the pro- 
tostellar environment, in particular by ultraviolet (UV) 
photons. However, this assumption has been challenged on 
the basis of a quantitative set of laboratory experiments 
studying the efficiency of OCN - photoprod uction, and the 
possibilities of forming OCN - the rmally I Raunier et alJ 
\200ii Ivan Broekhuizen et ai1l2004 - To assist the reader 
in following the subsequent discussions and conclusions of 
this paper, these two OCN" formation mechanisms, and 
the regimes in which they apply, are summarised below. 

In the laboratory, photochemical processes lead to 
the formation of OCN" at abundances of (at most) 
1.9% with respect to H^O-ice, starting from an in- 
terstellar ice analogue cont aining H2O, CO and NH3 
l|van Broekhuizen et alj[2004[) . The initial amount of NH3 
in these experiments ranges from 8 to 40% with respect to 
H2O ice. The abundance of interstellar NH3 is uncertain, 
but i s observed to be <5% toward W33A IjTaban et alJ 
120031) . Including the constraint that the abundance of 
NH3 ice after photolysis cannot be more than 5% low- 
ers the maximum OCN - photoproduction yield to 1.2% 
with respect to H?Q ice (s ee, for example, Figure 6 of 
Ivan Broekhuizen et alJl2004h . 

From these experiments, the maximum interstellar 
OCN - photoproduction yield was constrained to an 
abundance of 1.2% with respect to E^O-ice by as- 
suming a maximum abundance for solid NH3 of 5 % 
with respect to H ^O-ice, as observed toward W33A 
ijTaban et alJl20o3> . If the cosmic ray induced UV-fi eld 
(1.4 xlO 3 photons cm -2 s -1 . IPrasad fc Tarafdarl Il983l) is 
assumed to be the only source of UV-photons, this maxi- 
mum 1.2% OCN - abundance would be reached after a 
UV-fluencc equivalent to a molecular cloud lifetime of 
4xl0 8 years. This is long compared to the assumed age of 
the YSO's studied in this paper (~10 5 yr), even if a long 
pre-stellar phase of ~ 10 7 yr is included. 

Alternatively, laboratory studies show that OCN - 
can be formed via thermal heating of ices containing 
HNCO in the prese nce of NH3 or H?Q with an ef- 
ficien c y of 15-100% fPernvket all Il998l iRaunier et all 
2003; van Broekhuizen et alJ 120041) . Laboratory experi- 
ments and theoretical calculations show that this solva- 
tion process is even efficient at 15 K provid ed that HNCO 
is su f ficiently diluted in N H3 or H 2 0-ice ijRaunier et alJ 
l2003t IPark fc WoonlEool . for the proton transfer to oc- 
cur. 
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4. Decomposition of the XCN profile 



Our analysis of the XCN band, containing two components 
(see Sect. 12.2(1 . refines the results from Paper I where a 
single component fit was used. The compone nt centred at 
~ 217 5 cm" 1 has been discussed in detail bv lFraser et alJ 
( 2005) in terms of gas-surface interactions of CO with in- 
terstellar grains. Thus, in the remainder of this paper the 
discussion is focused on the carrier of the component of 
the XCN band centred at ~2166cm -1 , and simply uses 
a two component fit to correct for any contribution from 
the component centred at ~ 2175 cm . 

The distribution of the ^2166 cm" 1 component of the 
XCN band shows a broad spread of band centres, most 
probably related to source-to-source differences in the lo- 
cal ice-environment of the carrier of this component in the 
various lines-of-sight. A good 'prototype' for this compo- 
nent is the XCN band of W33 A, as it is also one of the 
most red-centred and strongest of the objects studied. It 
matches a Gaussian profile, peaking at 2165.7 ±0.1 cm -1 
(FWHM = 26 ± lcm" 1 ), that is well reproduced by lab- 
oratory spectra of OCN" (see Fig. OJ. An initial analysis 
showed that this Gaussian profile also reproduces the red 
wing and component of the XCN band of all other YSO's 
in this study. Consequently, the XCN band was analysed 
using a si milar proced u re to that employed for the 6.8 /xm 
band by iKeane et al ] l|200lh . whereby the 6.8 /xm band 
was also decomposed by using the most extreme source 
in the sample, MonR2:IRS3, to fit all lines-of-sight to- 
gether with a 'rest-spectrum' obtained after the initial fit- 
ting of the other sources. Here, the Gaussian fit to the 
W33 A XCN band is adopted (henceforth referred to as the 
2165.7-component), and not the actual YSO spectrum, to 
avoid any contribution from the CO-ice band to the fit. 

The component underlying the distribution of XCN 
bands centred at ~2175cm _1 in Fig. [2] was determined 
empirically, in effect by fitting the residual of the XCN 
band after spectral subtraction of the contribution of the 
2165.7-component. This resulted in a Gaussian, centred 
at 2175.4 ±3 cm" 1 (FWHM = 15 ±6 cm" 1 ). This feature 
will be referred to as the 2175.4-componcnt. Adopting a 
different Gaussian for the 2165.7-component, i.e. one cen- 
tred at 2166.9 or 2168.7cm -1 , does not significantly affect 
these fit results, which indicates that the range of environ- 
ments observed for the 2165.7-component do not introduce 
an error on the fit. 

The fitting routine is based on the same IDL algo- 
rithms developed in Paper I. All spectra were fitted us- 
ing the same 2165.7- and 2175.4-components, keeping the 
centre positions and FWHM constant. The optical depths 
of the two components are the only free parameters (see 
Tablc|21for a summary). In addition, a full fit of the CO-ice 
feature was run in parallel to correct for possible overlap 
between the wing of the CO-ice feature and the XCN band 
(see Paper I and Fig. 4b for details on the CO-ice fit). 



(a) 



HAFGL 70005 




(b) 




RAF GL 7009S 



10 

N(HyO) / 10 18 cm- 3 



14 

7 12 
" 10 



(C) 



RAFGL 70095 




N(H a O) / 10 18 cm- 2 

Fig. 5. (a) The optical depth of the 2165.7-component, 
(b) the optical depth of the 2175.4-component, and (c) 
the total integrated area of the XCN band, as determined 
by Eq.l, are shown versus the H2O ice column density, 
A r (H20). High-mass objects are indicated by diamonds, 
all others, dominated by low-mass YSO's, by stars. 3 a 
uncertainties are given, or upper limits are indicated by 
arrows. 



5. Results 

5.1. Results of the fitting procedure 

Fig. 2] shows the results of fitting the XCN band in five 
of the YSO's in our sample. These are presented from top 
to bottom in order of decreasing wavelength of the XCN 
band centre position. The XCN band of W33 A (Fig. Hi) 
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shows the excellent fit of the 2165.7-component with < 1 % 
contribution by the 2175.4-componcnt. In Fig. 0Jd, al- 
ready some contribution to this band from the 2175.4- 
component is apparent. From Fig. 0Jd to e, the 2175.4- 
component contribution increases, blue-shifting the XCN 
band, until it dominates the XCN band toward Elias32 
and IRS 63. Table ^ summarises the derived optical 
depths. The % 2 of the fit (not shown) generally varies be- 
tween 0.3 and 10. The additional three-component fit of 
the solid CO feature (Paper I) is shown by the three grey 
dash-dotted curves in Fig. 0b only, but was applied to all 
fits. 

5.2. Correlations with H2O 

In Fig. Et, the optical depth of the 2165.7-component 
is plotted against the H2O column density, A^^O). 
For all YSO's studied here, the derived optical depth 
of the 2165.7-component is <0.16, except for HH46, 
RAFGL 7009S and W33A. No relationship with N(R 2 0) 
is apparent. In fact, it seems that the optical depth value 
of the 2165.7-component toward high-mass YSO's is in 
most cases larger than toward lower mass YSO's, where 
the optical depth is typically < 0.09. 

Fig. Et> presents the optical depth of the 2175.4- 
component plotted against A^H^O). The optical depth 
of the 2175.4-component is <0.14 for all sources in this 
study and is not correlated to the H 2 col umn density. 
An in terpretation of this feature is given by iFraser et alJ 
ll2005h . 

Comparing Figs. Et and Et> shows that the 2165.7- 
component dominates the XCN band toward most high- 
mass YSO's, whereas the 2175.4-component contributes 
most to the XCN band in all other sources. Checks of a 
possible relationship between the optical depth of the two 
components (not shown) proved negative, indicating that 
the components do not trace the same feature. 

In Fig. Et the total integrated area of the XCN band 
is plotted against Af(H 2 0). The total integrated area of 
the XCN band is derived from the sum of the integrated 
areas of the individual components 



r (XCN) 



dl>= I T2165.7 dv 



T2175.4 dv 



(1) 



The XCN band is evenly distributed toward high- and low- 
mass YSO's and its total area tentatively increases with 
increasing A^(H20). 

5.3. Correlations with CO 

The CO-ice band, observed at around 2140 cm" 1 toward 
all lines-of-sight here, has been analysed in detail in 
Paper I, using a phenomenological decomposition into 3 
components to fit its band profile. The most red-centred 
component (rc) was attributed to CO in a hydrogen- 
bonded environment, although this cannot fully explain 
its extended red wing. It was suggested that the rc may 
evolve from pure CO-ice when thermal annealing of the ice 



n. 0.3 r 
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(c) 



RAFGL 7009S 



W33A x 4y HH 46 




RCi'A IRS7B 



1.0 0.8 0.4 0.6 0.8 1.0 1.2 1.4 

Optical depth rc of CO — ice 

Fig. 6. (a) The optical depth of the 2165.7-component, 
(b) the optical depth of the 2175.4-component, and (c) 
the total integrated area of the XCN band, as determined 
by Eq. 1 , are presented versus the optical depth of the red 
component (rc) of CO-ice. High-mass objects are indicated 
by diamonds, all others, dominated by low-mass YSO's, 
by stars. 3 a uncertainties are given, or upper limits are 
indicated by arrows. The shaded area, shown in each plot, 
marks the correlation between the optical depth of the 
XCN band and the red component (rc) of CO-ice as pro- 
posed by Paper I. 



induces CO to migrate into, and get trapped in, porous 
H2 0-ice. As such, this component may tr ace the ther- 
mal history of the interstellar ice mantle l(Tielens et alJ 
Il99ll Paper I). Hence, the relationship between the opti- 
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cal depth of the rc of CO-ice and the XCN band may pro- 
vide information on the thermal history of the XCN band 
as well. Paper I found a tentative correlation between this 
rc and the optical depth of the XCN band and suggested 
that this was actually due to a component of the XCN 
band, centred between 2170-2180 cm -1 . 

Fig. HJi presents the optical depth of the 2165.7- 
component versus the optical depth of the rc of the CO-ice 
band, showing that this component and the rc of CO-ice 
are not correlated. The shaded area marks the correlation 
suggested by Paper I. However, those sources that show a 
strong 2165.7-component also tend to have a strong CO 
rc, but not the other way around. 

In Fig. m> the optical depth of the 2175.4-component 
is plotted against the optical depth of the rc of the 
CO-ice band. In contrast to the 2165.7-component, the 
2175.4-component matches the correlation indicated by 
the shaded area much better, particularly at rc optical 
depths of below 0.6, corroborating the idea posed by Paper 
I. However, the YSO's with deep rc optical depths (and 
deep 2165.7-component optical depths) fall outside the 
correlation proposed by Paper I. 

For completion, Fig. HO: illustrates the relation between 
the integrated area of the XCN band and the optical depth 
of the rc of CO-ice, showing that most low-mass YSO's 
match the correlation suggested by Paper I (shaded area) 
and reflecting that in these objects the 2175.4-component 
dominates the XCN band (see Sect. I5.2|) . Conversely, ap- 
proximately half of the XCN bands observed toward high- 
mass YSO's clearly exhibit larger integrated areas than 
would be expected from that correlation. 

6. Discussion 

6.1. The OCA/ - abundance 

Based on the discussion in Sect. 13.11 only the 2165.7- 
component of the astronomical spectra can be confi- 
dently assigned to OCN - . Thus, the OCN - column den- 
sity, 7V(OCN - ), and its abundance with respect to H2O, 
[OCN - ], are derived from the optical depth of the 2165.7- 
component via: 



[OCN - ] = 



jy(ocN - ; 

W(H 2 0) 



j T2165.7 dv x Aq^ 



OCN- 



iV(H 2 0) 



(2) 



where the FWHM of T2165.7 i s fixed at 26 cm -1 and Anr.N- 
is 1.3xl0 -16 cmmolec. -1 ( van Broekhuizen et alJ 12004. 
see Sect. 13.21 for a discussion on the value of Aqcn-)- 
The inferred iV(OCN - ) and [OCN - ] are summarised in 
Table H 

In Fig. [7| the OCN - abundance with respect to 
iV(H 2 0), [OCN - ], is plotted against iV(H 2 0). All sources 
studied here show abundances of <1.0%, except for 
W33 A. Excluding the high-mass stars, gives abundances 
of < 0.85 %. The high signal-to-noise of the spectra allows 
the detection of the XCN band down to optical depths of 
0.01. As a result, very strict upper limits have been put 
on the OCN - abundance along some lines of sight. 



1.00 



0.10 



0.01 




W33A 



RAFGL 7009S 



N(H 2 0) /'10 18 cm- 2 

Fig. 7. The OCN - abundance, [OCN - ], plotted with re- 
spect to the H 2 column density, _/V(H 2 0). High-mass ob- 
jects are indicated by diamonds, all others, dominated by 
low-mass YSO's, by stars. 3<r uncertainties are given, or 
upper limits are indicated by arrows. The black horizontal 
line marks the 1.2% maximum abundance of OCN - that 
may be formed photochemically under interstellar condi- 
tions, as determined from laboratory studies. 

6.2. OCA/ - abundance variations 

The sample of XCN observations presented here is large 
enough to obtain statistically significant information on 
the distribution of abundances and to compare between 
different classes of cloud environments, such as low- and 
high-mass star-forming regions. As seen in Fig. [7| and 
Table0] the abundance of OCN - is observed to vary by at 
least a factor of 10-20 among the lines of sight in the sam- 
ple. This observed range of OCN" abundances is larger 
than that of CH3OH and second only to that of ice species 
whose abundances are sensitive to low-temperature freeze- 
out and evaporation, such as CO. 

A particularly interesting point is that some of the ob- 
served lines of sight toward low-mass sources are located 
in close proximity to each other in the plane of the sky 
(100-10000 AU). Such objects provide a direct probe of 
the physical size of regions with enhanced (or reduced) 
OCN - abundances. In particular, several close binaries in 
the sample probe variations in abundances on scales of a 
few hundred AU, i.e. at the scale of circumstellar disks. 
Some of the best upper limits on the OCN - abundance of 
< 0.1 % are observed toward IRS 43 and IRS 44 in the p 
Ophiuchus star-forming cloud. These lines of sight are lo- 
cated approximately 4' or 30 000AU from IRS 51, which 
has a measured OCN - abundance of 0.8%. This shows 
that OCN - is formed in localised regions even within the 
same low-mass star-forming cloud complex. The individ- 
ual components of the close binary source EC 90A+B, 
which were observed separately, appear to have OCN - 
abundances varying by a factor of 6. Since the separation 
between the two components is only 1'.'6, corresponding to 
approximately 400 AU, this is an indication that processes 
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Table 4. Results of the XCN band fitting procedure 



Source a 


t(2165.7) 6 


t(2175.4) 6 


JV(OCN~) c 


[OCN-] c 


t(H 2 0) " 




2165.7 cm -1 


2175.4 cm -1 


X 10 16 molcc. cm -2 


% 




p Ophiuchus 


IRS 42 


<0.023 


0.020±0.010 


<0.46 


<0.25 


1.17 


IRS43 


<0.02 


0.072±0.030 


<0.40 


<0.18 


1.40 


IRS44 


<0.02 


<0.024 


<0.40 


<0.13 


1.95 


IRS46 


<0.039 


<0.035 


<0.78 


<0.68 


0.74 


IRS51 


0.073±0.013 


0.037±0.015 


1.46±0.26 


0.76±0.12 


1.23 


IRS63 


<0.015 


0.075±0.016 


<0.3 


<0.17 


1.12 


WL12 




0.062±0.032 






1.95 


CRBR 2422.8 


<0.029 


0.120±0.045 


<0.58 


<0.19 


1.90 


Eliaa 32 


<0.030 


0.074±0.036 


<0.60 


<0.36 


1.04 


VSSG17 




<0.095 






0.87 


R.NO 91 


0.087±0.010 


0.066±0.017 


1.74±0.2 


0.40±0.09 


2.88 



Serpens 



EC90 A 


0.010±0.006 


0.039±0.008 


0.2±0.12 


0.15±0.05 


0.86 


EC90 B 


0.056±0.012 


0.030±0.007 


1.12±0.24 


0.86±0.12 


0.84 


EC82 


<0.02 


<0.026 


<0.40 


<0.86 


0.30 


CK2 


<0.027 


<0.089 


<0.54 


<0.19 


1.86±0.19 e 


SVS 4-9 


<0.08 


<0.04 


<1.60 


<0.64 


1.6±0.l'' 


SVS 4-5 


<0.073 


<0.010 


<1.46 


<0.27 


3.50±0.30 / 



TPSC 78 


0.155±0.003 


<0.025 


3.1±0.06 


1.12±0.14 


1.77 


TPSC 1 


0.158±0.032 


<0.042 


3.16±0.64 


0.97±0.13 


2.08 


Rcipurth 50 


0.100±0.005 


0.049±0.012 


2.0±0.1 


0.38±0.08 


3.42 


Corona Australis 


HH100 IRS 


0.04±0.002 


0.098±0.003 


0.8±0.04 


0.32±0.08 


1.6 


HH46 


0.245±0.019 


0.095±0.032 


4.9±0.38 


0.59±0.01 


5.3±3.0 9 


RCrA IRS7A 


<0.029 


<0.074 


<0.58 


<0.11 


3.5 


RCrA IRS7B 


<0.035 


<0.050 


<0.70 


<0.13 


3.5 


RCrA IRS5A 


0.040±0.025 


0.123±0.036 


0.80±0.5 


0.20±0.06 


2.60 


RCrA IRS5B 


<0.022 


<0.054 


<0.44 


<0.12 


2.42 


Chamaeleon 


ChaIRN 


0.088±0.069 


<0.055 


1.76±1.38 


0.75±0.12 


1.5 


ChaIRS 6A 


<0.05 


<0.108 


<1.0 


<0.22 


2.90 


Vela 


LLN17 


0.104±0.011 


0.048±0.014 


2.08±0.22 


0.62±0.11 


2.14 


LLN20 


0.087±0.040 


<0.024 


1.74±0.8 


0.41±0.09 


2.75 


LLN33 


0.130±0.032 


<0.051 


2.6±0.64 


0.62±0.11 


2.69 


LLN39 


<0.02 


<0.02 


<0.40 


<1.3 


0.20 


LLN47 




0.050±0.033 






1.15 


Taurus 


LDN 1489 IRS 


<0.014 


0.070±0.016 


<0.28 


<0.06 


3.18 


Additional sources ft 


RAFGL 2136 


0.123±0.011 


<0.051 


2.46±0.22 


0.47±0.09 


3.33' 


NGC7538 IRS1 




<0.019 






1.47 


RAFGL 7009S 


0.408±0.021 


0.104±0.026 


8.16±0.42 


0.71±0.11 


7.33 J ' 


W33A 


1.204±0.030 


<0.004 


24.08±0.60 


1.93±0.01 


8.00±0.75 1 ' 


RAFGL 989 


0.026±0.013 


0.019±0.017 


0.52±0.26 


0.23±0.06 


1.45 



a Source spectra taken from Paper I, unless otherwise stated. 3<r uncertainties or upper limits are given. c Column densities 
and abundances derived from the 2165.7-component using Eq.2. d H2O optical depths from Dartois et al. private communi- 
cation , u nless otherwise stated, and d erived according to Se ct. 12. II Uncertainties ar e between 10 to 20%. e lEiroa Hoo^aprJ 
1 19891) f iPontoppidan et al] <2004t) 9 iBooeert et aiTj2004) h ISO-data archive. i iGerakines et alJ lll999l) j iDartois et alJ 
1 19991) fc iGibb et all J2000|) 



occurring on the size scale of circumstellar disks may be 
required to form OCN - in large quantities. Note, however, 
the spectrum observed toward EC 90B is complicated 
by complex emission and absorption from gaseous CO, 
so that the observed difference in OCN - abundances for 
this binary requires further investigation. The two com- 
ponents in another close binary, RCrA 1RS5, also show a 
tentative difference in OCN - abundance of a factor of 2. 
Clearly, only slightly more sensitive observations of em- 
bedded close binaries will be able to confirm the obser- 
vation of a varying OCN - abundance on scales of a few 
hundred AU. 

6.3. OCA/ - formation toward low-mass YSOs 

The observed abundances and their variations for our large 
sample of sources put important constraints on the forma- 
tion mechanism of OCN - . All inferred OCN - abundances 



presented here allow quantitatively for a photochemi- 
cal formation mechanism (see Sect. Only W33A, 
an embedded O-typ e star with a luminosity of 10 5 L© 
l|Mueller et all2 002) with an OCN - abundance of ~ 2% is 
at the edge of the range, as are the Galactic center source 
GC:IRS19 and the external galaxy NGC 4945 presented 
elsewhere (see van Broekhuizen et al. 2004). If we define 
high- and intermediate mass YSO's as those more lumi- 
nous than 50 i©, it is found that OCN" abundances of 
order 0.85-1% are observed toward high- intermediate- 
and low-mass stars in our sample, but strict upper limit 
abundances of < 0.2 % are observed only toward low- 
mass stars. This does indicate that some differences ex- 
ist in the efficiency of OCN - production between high- 
and low-mass star-forming regions, and may suggest that 
UV photons play a role. However, comparison with lab- 
oratory studies shows that the photon dose required to 
produce the highest OCN - abundances (of ~ 0.85 %) de- 
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tected here toward low-mass YSO's are ~ 40 x larger than 
can be accounted for by the cosmic ray induced photons 
that any ice might receive within the typical lifetime of the 
pre- and protostellar phases (see Sect. l3.3|l . Adopting the 
somewhat higher cosmic-ra y induced photon fluxes esti- 
mated bv lShen et alJ l)2004|) would lower these timescales 
by an order of magnitude, and they could be further de- 
creased by a factor of 2 if the initial NH3 abundance were 
as high as 30%, but it remains difficult to fully reconcile 
the ages. 

Another, perhaps more serious problem with the UV 
photoproduction scenario is that little variation in the 
OCN - abundance from source to source is expected if 
the cosmic-ray induced photon flux is assumed to be the 
same everywhere. Thus, local sources of radiation — either 
UV from the young star or star-disk boundary — need to 
be invoked to explain the observed small-scale variations. 
A well-known concern is whether any of this UV radia- 
tion is able to reach the bulk of the ices because of the 
large extinctions involved. Spherically-symmetric models 
of YSO envelopes have been constructed, taking the vari- 
ation of the density, dust te mperature and UV rad iation 
with depth into account (e.g. IjQrgensen et"aLll2002|) . and 
it has been shown that any UV from the st ar only reaches 
the in ner few % of the envelope mass (e.g., IStauber et all 
2004). Thus, special geometries are required to allow the 
UV to escape and reach the bulk of the ice, such as out- 
flow cavities or flaring turbulent disks in which the ices 
are cycled from the shielded midplane to the upper lay- 
ers. One alternative option is UV radiation created by 
secondary electrons resulting from X-rays from the young 
star. Model calculations show that the X-rays, and thus 
the UV radiation, penetrate much deeper into the enve- 
lope than the UV from the star itself so that it can reach 
a larger fraction of the ices l)Stauber et aLlEoOfil) . It would 
be worthwhile to couple such models with an ice chemistry 
to check the OCN - production quantitatively, but this is 
beyond the scope of this paper. 

An alternative scenario for OCN - formation is through 
thermal heating of HNCO in the presence of NH3- or 
fi20-ices. Assuming a 15-100% efficiency of this mech- 
anism as determined under laboratory conditions (see 
Sect. 13. the OCN - abundances toward low-mass YSO's 
presented here imply that 0.85-5% HNCO should be 
present in the solid state prior to the formation of 
OCN - to explain the highest observed abundances of 
-0.85%. To date, HNCO has been detected in inter- 
stellar gas associated with star-forming regions, with 
column densities in the range from 0.3-7.5 x!0 14 cm~ 2 
llWang et alJl2004 IZinchenko et allEoOot iKuan fc Snvderl 
ll996H,Iackson et al.lll984Uvan Dishoeck et alJll995h . This 
is a factor of 100 less than what is expected from 
the OCN - column densities derived here, assuming all 
the OCN - evaporates as HNCO. In the solid state, 
only HNCO upper limits, of the order of 0.5 to 0.7% 
with respect to H 2 0-ice, are dete rmined toward W33A, 
NGC7538 IRS9 and AFGL2136 l|van Broekhuizen et alJ 
l2004|) . Toward the low-mass YSO Elias29, this is < 0.3% 



(derived from lBoogert et allEoO^ . The ease with which 
HNCO is deprotonated in a H20-rich environment implies 
that HNCO would be difficult to detect toward regions 
where the ice is warmer than 15 K, although its presence 
may be observable by means of its 4.42 feature toward 
the coldest regions where the ice temperature is close to 
10 K. On the basis of tqhn- derived here, and t he band 
strength of HNCO l|van Broekhuizen et "aT1l2004[K t H nco 
is expected to range from 0.001 to 0.08, a factor of 5 less 
than observed for OCN - . 

Chemical models predict that HNCO is formed on in- 
terstellar grains via grain-surface reactions with an abun- 
dance upper limit of < 1-3 % with respect to H2 0-ice 
l|Hasegawa fc Herbstlll993l iKeane et al.ll200lh . Given the 
observed upper limits and model abundances for HNCO, 
thermal processing also has difficulties to explain the high- 
est observed OCN - abundances, but it is a serious alterna- 
tive to UV-photoprocessing for the formation of at least a 
fraction of OCN" in interstellar ices. The observed OCN - 
abundance variations would in this scenario have to result 
from variations in the grain-surface chemistry rates form- 
ing the reactants that lead to OCN - , for example due to 
temperature or availability of atomic N. 

On the basis of its spectroscopy, it is not possible to 
draw a distinction between OCN - formation via proton- 
, electron-, UV- or thermally induced processes, neither 
from laboratory spectra fSect. l3~T)l . nor from astronomi- 
cal observations of the 2165.7 cm -1 OCN - band. 

7. Conclusion 

The XCN band has been observed toward 39 YSO's en- 
abling the first detailed study of the OCN - abundance 
in a large sample of low-mass YSO's. Statistical analysis 
of the band centre position distribution proves unequiv- 
ocally that the band contains at least two components, 
of which only one (centred at 2165.7cm -1 ) can be associ- 
ated with OCN - as its carrier, based on laboratory OCN - 
spectra. We conclude that in all cases the OCN" is em- 
bedded in a strongly hydrogen-bonded, and possibly ther- 
mally annealed, ice environment. A phenomenological de- 
composition was undertaken to fit the full XCN band pro- 
file toward each line of sight using two components, cen- 
tred at 2165.7cm- 1 (FWHM = 26cm -1 ) and 2175.4cm -1 
(FWHM = 15cm -1 ). OCN - abundances were derived 
from the 2165.7-component of this fit. 

Typically, OCN - is detected toward the low-mass 
YSO's in our sample at abundances <0.85% (with re- 
spect to the H20-ice column density) and < 1.0% toward 
high-mass YSO's, with the exception of W33 A. A large 
dynamic range of abundances is observed, varying by at 
least a factor of 10-20. Together with the variations in 
abundance observed between regions separated by only 
400 AU, this provides a strong indication that OCN - for- 
mation is a highly localised process. The inferred OCN - 
abundances allow quantitatively for a photochemical for- 
mation mechanism, but the observed large variations are 
more difficult to explain unless local sources of UV radia- 
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tion or X-rays are invoked which reach a large fraction of 
the ices. Alternatively, OCN - can be formed in the bulk 
of the ice from the solvation of HNCO, which may have 
formed via surface reactions. However, this too has diffi- 
culties. Further diagnostics are required to be more con- 
clusive about the chemistry leading to OCN - formation in 
interstellar environments. These include the construction 
of more detailed physical models of the sources presented 
here to derive the effects of varying grain temperatures 
and UV radiation fields on the different formation mech- 
anisms of OCN", slightly more sensitive observations to 
probe OCN - abundance variations on a 100 AU scale, and 
laboratory data on the formation of HNCO via surface re- 
actions at various temperatures. 
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